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Abstract 

The sources of the systematical uncertainties due to radiative correction proce- 
dure in measurements of the structure function with 3 He target are considered. 
Their numerical estimations are presented. The relative systematical uncertainty 
does not exceed 5%. 



The measurements of the spin-dependent structure functions performed in the last 
years (SMC, E154/155) have the tendency to decrease essentially the statistical uncer- 
tainties comparing with previous ones (EMC, E142/143). Thus the one of main problem 
is to decrease the magnitude of the systematical uncertainties or define more correct way 
of their calculation. In general systematical uncertainties come from measured quanti- 
ties, namely: values of the beam and target polarizations, efficiencies of the coordinate 
detectors etc. There are also some contributions in total systematic uncertainty due to 
unmeasured quantities such as radiative corrections (RC) to be applied to extract the 
one-photon exchange cross section from the measured one. 

So called radiative events, which originate from loop diagrams and from processes 
with the emission of additional real photons, cannot completely be removed by exper- 
imental methods and so they have to be calculated theoretically and subtracted from 
measured cross sections. The calculation of the radiative corrections requires knowledge 
of spin-independent and spin-dependent sturcture functions both in region measured in 
the considered experiment and beyond it. The choice of different parameterizations of the 
sturcture functions, elastic and quasielastic formfactors, neglecting of electroweak and 
higher oreder effects and simplifications in RC procedure leads to uncertainties in calcula- 
tions of RC. The approach of calculation of these systematical uncertainties is presented 
in this report. We consider the case of the 3 He target and kinematical region close to 
HERMES IJ and SLAC 0j experiments. To get the numerical estimations of the system- 
atical uncertainties due to radiative correction procedure we use the special program || 
and radiative correction calculation code POLRAD 2.0 (|, |5|. 

The extraction of the spin asymmetry A\(x, Q 2 ) in the experiments measuring spin- 
dependent structure functions is based on the following formula 

1 jv Ti - iV TT 

Mx ' Q) = <Df d > PTW' (1) 

where <Df d > is mean value of product of depolarization and dilution factors (see |J for 
details), N*' and iV^J- are the number of events for parallel and antiparallel spin target 
configurations. 



In order to take into account the radiative effects these numbers have to be calculated 
as the weighted sum 

N n,n = ^ w(x,Q 2 ). (2) 
tut 

Here the weight (w) is calculated as a ratio of Born and observed cross sections 

w = ■ (3) 

O"0 + &RC 

We will refer to this procedure as 'exact'. 

The exact procedure requires the RC calculation for each event, so in practice other 
scheme is used. We will refer to this procedure as 'standard'. In this case the RC is 
applied to asymmetry averaged at x bins. 

The radiative correction AA\ to the measured asymmetry is defined as: 

Al meas = A X + AA 1 (4) 

and can be written in terms of spin- independent (a") and spin-dependent (a p ) parts of 
DIS cross section. 

AA _ ^KM + g + <d - gg(gi)(g& + g + gg) (5) 

where 5 V = &p/(Jq = o""/ctq. The polarized parts of Born cross section and inelastic 
radiative tail depend on g\. 

The radiative correction procedure is performed as follows. 

1. The measured asymmetry A r y i (i = 1, N x , where N x is number of x bins) is fitted 
by a functionF] with taking into account statistical uncertainties of A™. 

2. The constructed fit is used for calculation of of n (gi) and Oo(gi). 

3. The extracted asymmetry is calculated for each kinematical bin as follows 

AT! = Aft - AA, (6) 

There are three important sources for uncertainties coming from radiative correction 
procedure: a) using of simplified (standard) scheme instead of the exact one ; b) using 
of models and data for structure functions; c) physical effects which are neglected in the 
standard scheme. 

To calculate the systematical uncertainties of the types a), b) and c) the following 
scheme is used. The Monte-Carlo kinematical events are generated according to random 
flat generator: In Q 2 for Q 2 over allowed Q 2 - region and flat for v over allowed v - region. 
After calculation of kinematic variables and applying of kinematic cuts close to acceptance 
of the HERMES experiment [IJ, the weights with DIS cross section u> l7 are obtained at 
given kinematical point. The model for spin asymmetry and generated Born asymmetry 
are plotted on fig.pi. 

For each event the recalculation of the weight is performed with radiative correction 
factor: 

w° bs = — , (7) 
w 



1 Th.e example such a function for 3 He target can be found in Appendix of ref 
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where weight w defined from eq. @ is calculated by code POLRAD 2.0. The measured 
(A° bs ) and extracted within exact scheme (Af) asymmetries are obtained using formulas 
(H) and (H) with the weights w obs and w 11 . 

The application of the standard scheme to A° bs gives the extracted asymmetry Af 
within standard scheme. To estimate the relative systematical uncertainty of type a), 
we calculate the difference between extracted asymmetries within the exact (Af) and 
standard {Af ) schemes for each kinematical bin: 



| yje^ A | 

^ a A ex 

The relative systematical uncertainties of types b) and c) are estimated as follows 



Af 



(9) 



The asymmetry Af is calculated using the weight w ex = w w. The weight w is also cal- 
culated according to formula @, but a obs is calculated with different models for structure 
functions or with taking into account electroweak and high order effects. 
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Figure 1: The generated asymmetry along used fit and relative systematical uncertainties 
versus x in percents due to simplification procedure. It should be noted that peak for 
x ~ 0.4 is artificial and originates from the fact that neutron spin asymmetry is close to 
zero in this kinematical point. 

The quantity e a can be non-zero only due to order of averaging within standard and 
exact schemes for the same models for structure functions. In case of 'standard' scheme, 
the results are calculated for averaged kinematical point at each bin. But this mean value 
can be shifted due to very different dependence of Born and radiative corrected cross 
sections on kinematical variables. The exact scheme is free from this shortcoming. In 
this case the RC is taken into account before averaging using calculation of the weight (|3]) 
for each event. The relative systematical shift obtained with standard and exact schemes 
does not exceed 2% (fig.|l|b). 

The quantity (Tq in eq.(||D depends on F2 in the given kinematical point Q 2 and x, but 
o~i n (see eq. (|J)) requires knowledge of the structure function F 2 in the wide region of 
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varying of kinematical variables x and Q 2 . So the fit of F 2 used for RC calculation has to 
describe adequately both resonance and DIS region as well as to have correct asymptotics 
behaviour for Q 2 — * and W 2 — ► (M + m n ) 2 . Such fit can be constructed on the basis of 
NMC parameterization of F 2 for protons and deuterons 

F?\x } Q 2 ) = ^F 2 H (x,Q 2 )+F 2 D (x,Q 2 ))j (10) 

F 2 and F 2 are taken from fit described in 0. The fit takes into account the contribution 
of A resonance, has correct behaviour on boarders and describes DIS data for x > 0.01. 

Another possibility is to change the model for F 2 on standard POLRAD fit ||, which 
includes Brasse parameterization of the three resonances (instead of one resonance as 
in standard fit (0)), Stein fit for small Q 2 region and 15 parameters NMC fit for DIS 
region. The final results shows that unpolarized structure function can change asymmetry 
by approximately 1%. 

Both polarization and unpolarization contribution of quasielastic radiative tails depend 
on quasielastic response functions Tl (see [|J for details). The electric and magnetic 
formfactors for proton and neutron fall as Q 8 for high Q 2 . So the only region of small 
Q 2 ~ M 2 is important. In this region the nucleon formfactors are known with good 
accuracy, and their variation does not lead to systematical uncertainty. To the contrary 
the models for electric S#, magnetic Sm and mixed Sem suppression factors differ in this 
region. The code POLRAD 2.0 exploits Y-scaling hypothesis || |9] Sm = Se = Sem = 
F(uq) and scaling function F(vq) are calculated in Fermi gas model |T0|, [TT|] . Alternatively 
the suppression factors can be also calculated within the sum rule approach ||12|| . The 
difference is important for small values of x, which correspond to high y, where it can 
reach 1.5%. 

The elastic structure functions are calculated as quadratic combinations of electric 



and magnetic formfactors. A simple Schiff's model with gaussian wave function [T3j is 
used as an alternative model. The results are similar to quasielastic case. 

Resonance region gives a large contribution to RC for spin- independent DIS. The 
contribution to resonance region in g% can be also important. Unfortunately, there are 
no enough experimental data or satisfactory models for g™(x, Q 2 ) in resonance region 
(W 2 < 4). It is the main reason why scaling behaviour of spin asymmetry is extrapolated 
into resonance region under POLRAD consideration. For alternative approach we used 
simple model for structure function g™(x) in region of A(1232) resonance constructed 
on the basis of two assumptions, namely: the iy 2 -dependence has the Gaussian form 
with height and width which can be roughly estimated from recent SLAC data UM, 



the Q 2 -dependence is defined by resonance contribution to Drell-Hearn-Gerasimov sum 
rule given in ref.[[HJ. The relative systematical uncertainty due to resonance region is 
important (~1.5%) for small x-region. 

The kinematical coefficient front of structure function g 2 is small enough for both at 
the Born level and for RC cross section. So, the contribution of the structure function is 
neglected normally. To study influence of such an approach the model of the Wandzura- 
Wilczek is applied ([T1J . The small effect is obtained. The exception is last two bins where 



e b ~ 1%. 

The electroweak effects are not included in standard radiative correction procedure, 
because for current polarized experiments Q 2 ~ lOGeV 2 <C M 2 (M z is the Z-boson mass) 
and hence their contribution is small, but it has to be added to systematical error. Such 
a systematical uncertainty is estimated at the born level using code POLRAD 2.0. The 
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electroweak correction cannot be calculated by model independent way, that is why the 
quark parton model was used and GRV- , GRSV-parameterizations JT7|, [L8| for spin- 
independent and spin-dependent partonic distributions were applied. The correction is 
important (e c ~ 1.5%) for high x. 

In standard consideration the effects of higher order is estimated by simple exponen- 

19]. The POLRAD 2.0 gives also a possibility to 

The radiation 
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tiation procedure of soft photons 

obtain the a 2 order correction within structure function approach 
is considered to be collinear. There are three possibilities: initial, final state radiation 
and contribution of the Compton process. The first and second are important for inelastic 
radiative tail. The last one is extremely important for elastic and quasielastic radiative 
tails. The relative systematical uncertainties due to the higher order effect can exceed 2% 
at low x-region. 

The systematical uncertainties considered in previous sections can be gathered to- 
gether to obtain the total relative contribution due to RC to systematical uncertainties. 
The following systematical uncertainties are considered to be independent 

• simplification of procedure 

• unpolarized structure function F 2 

• quasielastic structure function 

• elastic formfactors 

• resonance region of structure function gi 

• polarized structure function g 2 

• electroweak effects 

• higher order effects 
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Figure 2: The relative systematical uncertainties due to radiative corrections. 



The quadratic sum of the uncertainties mentioned above presents in fig.|2|a. For x > 0.3 
the total systematical error does not exceed 2%. However for the first x bins with x < 0.3 
and y ~ 0.85 the effect is larger and can reach 5%. 
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Note that the uncertainties of types c) can be rejected from total sum, if they are 
included in standard scheme. For higher order correction it can be done without additional 
assumption, but electroweak correction calculation requires usage of the quark-parton 
model. In this case systematical uncertainties due to RC come only from procedure 
simplification and uncertainties in structure functions. This result is presented on fig.^b. 
It does not exceed 4%. 
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